Mon. Not. R. Astron. Soc. 000, [TJ{8] (YYYY) 



Printed 1 February 2008 



(MN M£X style file v2.2) 



Radio emission from the Sy 1.5 galaxy NGC 5033 



M. A. Perez-Torres 1 * and A. Alberdi 1 

1 Instituto de Astrofisica de Andalucia, CSIC, Apdo. Correos 3004, E-18080 Granada, Spain 



Accepted 2007 May 4. Received 2007 May 2; in original form 2007 March 23 



ABSTRACT 

We present new continuum VLA observations of the nearby Sy 1.5 galaxy NGC 5033, 
made at 4.9 and 8.4 GHz on 8 April 2003. Combined with VLA archival observations 
at 1.4 and 4.9 GHz made on 7 August 1993, 29 August 1999, and 31 October 1999, we 
sample the galaxy radio emission at scales ranging from the nuclear regions ( < 100 pc) 
to the outer regions of the disk (~ 40 kpc). The high-resolution VLA images show a 
core-jet structure for the Sy 1.5 nucleus. While the core has a moderately steep non- 
thermal radio spectrum (S v cx v a \ a\'\ « —0.4), the inner kpc region shows a steeper 
spectrum (af | k —0.9). This latter spectrum is typical of galaxies where energy losses 
are high, indicating that the escape rate of cosmic ray electrons in NGC 5033 is low. 
The nucleus contributes little to the total 1.4 GHz radio power of NGC 5033 and, 
based on the radio to far-infrared (FIR) relation, it appears that the radio and far- 
infrared emission from NGC 5033 are dominated by a starburst that during the last 
10 Myr produced stars at a rate of 2.8 M yr _1 yielding a supernova (type Ib/c and 
II) rate of 0.045 yr^ 1 . This supernova rate corresponds to about 1 SN event every 
22 yr. Finally, from our deep 8.4 GHz VLA-D image, we suggest the existence of a 
radio spur in NGC 5033, which could have been due to a hot superbubble formed as 
a consequence of sequential supernova explosions occurring during the lifetime of a 
giant molecular cloud. 

Key words: Galaxies: individual: NGC 5033 - Galaxies: Seyfert, starburst - Radio 
continuum: general - Radiation mechanisms: non-thermal 



1 INTRODUCTION 

Radio emission from normal galaxies is dominated by syn- 
chrotron emission from relativistic cosmic ray electrons 
(CRE) at low frequencies ( < 30 GHz), while at high fre- 
quencies (v > 30 GHz) the dominant emitting mechanism 
is thermal free-free emission from ionized gas at tempera- 
tures of ~ 10 4 K. 

With the advent of the Infrared Astronomical Satellite 
(IRAS, Neugebauer et al. 1984), a sample of « 20, 000 galax- 
ies complete to 0.5 Jy in the 60 /im band was detected, the 
majority of which had not been previously cataloged. While 
the observed FIR luminosity from most of those galaxies 
can be explained by ongoing star formation and starburst 
activity, some of the most luminous infrared galaxies may 
host an active galactic nucleus (AGN). Nearly all of the 
radio emission at wavelengths longer than a few cm from 
such galaxies is synchrotron radiation from relativistic elec- 
trons and free-free emission from H II regions. In turn, this 
synchrotron radio emission is a direct probe of the recent 
massive starforming activity, since only stars more massive 
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than about 8M0 result in the Type Ib/c and II supernovae, 
which are thought to accelerate most of the relativistic elec- 
trons that are responsible for the observed synchrotron radio 
emission. This results in an extremely good correlation be- 
tween radio and far- in f rared emission among galaxies (e.g., 
Harwit fc Pacinil 1 19751 . ICondon et alj Il991al . ICondon et al.l 
1991bl h If an AGN is present, then the associated radio con- 
tinuum, in excess of the level expected from the star for- 
mation, may be detected and, by means of high-resolution 
observations, its contribution to the total radio emission sub- 
tracted. 

Synchrotron radio emission also carries relevant in- 
formation on the magnetic field and the particle feeding, 
or re-acceleration mechanisms present in galaxies. Since 
the synchrotron energy loss rate of relativistic electrons, 
dE/dt, varies with the particle energy squared, E 2 , the 
typical lifetime of a relativistic electron, subject to radia- 
tive synchrotron losses is t syll — E/ (dE/dt) cx E . Fur- 
ther, the critical frequency at which an electron radiates 
most of its synchrotron emission is v c cx E 2 , and hence 
dE/dt cx u c , and t Byn cx u c 1//2 . Therefore, regions with lit- 
tle ongoing relativistic particle re-supply, or re-acceleration 
have steeper synchrotron spectra, which may be revealed by 
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multi-frequency radio interferometry observations. The sites 
and length scales where the spectral evolution is happening 
can also be revealed by high-resolution radio interferometry 
observations, which tell us where CRE are accelerated. 

High-resolution radio observations also allow the de- 
tailed study of both thermal and non-therm al emission in 
nearby, normal galaxies (e.g. ICondonl Il992h . Those obser- 
vations are needed, as the spatially integrated nonthermal 
spectral index is in general a poor diagnostic for the type of 
propagation or the importance of energy losses. On the con- 
trary, spatially resolved radio data for the haloes of galaxies 
allow us to draw firm conclusions. For example, the steepen- 
ing of the spectrum away from the disk is an indication that 
synchrotron and inverse Compton losses are taking place 
during the prop agation of cosmic ray electrons in the halo 
l|Lisenfeld fc Voelkll2000D . 

NGC 5033 is a nearby (D=13 Mpc) spiral galaxy that 
exhibits strong CO emission o ut to at least R ~ 50'.', and 
contains a Seyfert 1.5 nucleus ijHo et al.l[l997l ^. NGC 5033 
also displays st rong far infrared lum inosity (£[8-iooom m] = 
1.20 x 1O 1O L ; ISanders et alj|2003h , and has been host to 
at least three supernovae in the last 60 yr (SN 1950C, 
SN 1985L, and SN 2001gd). Here, we present and dis- 
cuss continuum VLA-D observations of NGC 5033 at 1.4, 
4.9 and 8.4 GHz taken on April 8, 2003, complemented 
with archival VLA data taken on August 29, 1999 (1.4 
GHz, A-configuration) and October 31, 1999 (4.9 GHz, B- 
configuration) , respectively. The resolutions yielded by the 
different arrays and frequencies allow us to trace the syn- 
chrotron emission not only in the innermost regions of the 
galaxy, but also in the inner and outer disk, and the galactic 
halo. 

The paper is organized as follows: we report on the radio 
observations in Sect. [2] we present our results in Sect. [3] and 
discuss them in Sect. [4] We summarize our main conclusions 
in Sect. [5] We assume throughout the paper a distance of 
13 Mpc to the host galaxy o f SN 2001gd NG C 5033, based 
on its redshift (z = 0.002839; iFalco et aljl997l l and assumed 
values of Ho — 65 km s _1 Mpc -1 and qo~0- At the distance 
of NGC 5033 1" corresponds to a linear size of 64 pc. 



2 OBSERVATIONS 

We reanalysed archival VLA data on NGC 5033 ob- 
tained under program AU0079 on 1999 August 29 (1.4 
GHz, A-configuration) and 1999 October 31 (4.9 GHz, B- 
configuration) , respectively. These data sets provide very 
similar angular resolution, since the increase in angular res- 
olution from 1.4 to 4.9 GHz is essentially compensated by a 
decrease in the longest antenna separation, which is approxi- 
mately a factor of three between the A and B configurations 
of the VLA. In both observing epochs, the VLA recorded 
both senses of circular polarization. Each frequency band 
was split into two intermediate frequencies (IFs), of 50 MHz 
bandwidth each. 3C 286 was used to set the absolute VLA 
flux density scale, and J1310+323 was used as the phase 
calibrator at both frequencies. 

We also show unpublished archival 1.4 GHz continuum 
VLA data obtained under program AP270 on 1993 August 
7, when the VLA was in C-configuration. The VLA recorded 
only the left sense of circular polarization, and the observa- 



tions were made in spectral line mode, since the aim was to 
image the 21-cm hydrogen line emission of NGC 5033. The 
6.2-MHz IF at 1.4 GHz was split into 63 spectral channels, 
each of 97.66 kHz bandwidth. For the purposes of imag- 
ing the continuum radio emission at 1.4 GHz, we used the 
channel-0 data, which contains the inner 75% of the line 
data. 3C 286 was used to set the absolute VLA flux density 
scale, and the source 1323+321 was used as the phase cal- 
ibrator. While the total synthesized bandwidth was small, 
the large time on-source (about 5.5 hr) and excellent uv- 
coverage allowed us to obtain the deepest 21 cm image ever 
of the Sy 1.5 galaxy NGC 5033 (see top left panel in Fig. QJ. 

Finally, we show observations at 4.9 and 8.4 GHz of 
SN 2001gd and its host galaxy NGC 5033, made on 2003 
April 8 with the VLA in D-configuration, as part of our 
VLBI observin g program to monitor th e angular expansion 
of SN 2001gd (jperez-Torres et all [20051 ). The VLA recorded 
both senses of circular polarization, and each frequency band 
was split into two intermediate frequencies (IFs) , of 50 MHz 
bandwidth each. We used 3C 286 to set the absolute VLA 
flux density scale at both frequencies. We used J1310+323 
as the phase calibrator at 4.9 GHz, while J1317+34 was used 
as the phase calibrator at 8.4 GHz. 

For all data sets, we used standard calibration and hy- 
brid mapping techniques within the Astronomical Image 
Processing System (A IPS) to obtain the images shown in 
Figure [T] 



3 RESULTS 

Figure Q] and Table [T] summarize our results for the total 
intensity radio images of NGC 5033, obtained with the VLA. 
The flux density errors given for the VLA measurements 
in Table [T] represent one statistical standard deviation, and 
are a combination of the off-source rms in the image and 
a fractional error, e, included to account for the inaccuracy 
of VLA flux density calibration and possible deviations of 
the primary calibrator from an absolute flux density scale. 
The final errors, <r/, as listed in Table [T] are taken as cr 2 — 
(eSo) 2 + <t 2 where So is the measured flux density, ao is the 
off-source rms at a given frequency, and e=0.05 at 1.4 GHz, 
and 0.02 at 4.9 GHz and 8.4 GHz, respectively. We note 
that the total quoted flux density errors are dominated by 
the uncertainties in the VLA calibration (e). The positions 
of the peaks of brightness for all the VLA images coincide 
among them within the errors, and correspond to the nucleus 
of NGC 5033. 

The 1.4 GHz VLA-C image (panel a in Fig. 1) shows 
the whole disk of NGC 5033, with an angular size of about 
10' x 5' (~ 38kpc x 19 kpc). The low-surface brightness, 
extended emission clearly follows the optical emission from 
the spiral arms of the galaxy, and most of the local max- 
ima of both optical and radio emissio n are spatially coinci- 
dent with H II regions (Evans 1996). The total flux den- 
sity above the lowest drawn contour is about 355 mjy, 
of which ~58 mjy are within a circumnuclear region of 
30"x30"(fs 1.8 x 1.8 kpc 2 ). We also detected three strong 
sources embedded within the extended 1.4 GHz emission of 
the disk, which we dubbed Nl, SW1, and SW2 in panels (a) 
and (f), and whose nature is discussed in section [4] 

The extension of the 4.9 GHz VLA-D total intensity 
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Figure 1. (a) Contours of 1.4 GHz observations made on 7 August 1993 with the VLA-C, overlaid on a Digitized Sky Survey (DDS2) 
Blue plate of NGC 5033; (b,c) total intensity radio images of NGC 5033 obtained from 4.9 and 8.4 GHz VLA-D observations on 8 April 
2003; (d,e) 1.4 GHz VLA-A observations on 29 August 1999 and 4.9 GHz VLA-B observations on 31 October 1999; (f) blow-up of the 
1.4 GHz VLA-C image, around the SW1 and SW2 regions. VLA contours are drawn at (-3,3,3 \/2,9,...)x the off- source rms of each map, 
which are of 93 fijy beam" 1 (1.4 GHz VLA-C), 45/xJy beam" 1 (4.9 GHz VLA-D), 19 fijy beam" 1 (8.4 GHz VLA-D), and of 170 and 
70 /xjy beam -1 at 1.4 and 4.9 GHz (panels d and e), respectively. The peaks of brightness are of 31.8, 10.6, and 2.4 mjy beam -1 at 1.4, 
4.9, and 8.4 GHz (top panels), and of 5.7 and 3.4 mjy beam -1 at 1.5 and 4.9 GHz (panels d and e). The principal dimensions (major 
axis X minor axis, position angle) of the restoring beams are: 17'.'3 xl5'/4, 83°(a); 19'.'1 xl5'/9, 83°(b); lO'.'l x8'.'7, -61°(c); 2'.'3 xl'/2, 
78° (d); T/8 xl'.'2, 82° (e). 



emission of NGC 5033 (panel b) is about 230" x 130" (~ 
14.7 kpcx 8.3 kpc), significantly smaller than the 1.4 GHz ra- 
dio emitting region. Since the 1.4 GHz VLA-C observations 
have a very similar resolution to our 4.9 GHz VLA-D obser- 
vations, this finding indicates that the external regions of the 
disk are partially resolved at this frequency. The total flux 
density of the nuclear regions of NGC 5033 above the lowest 
drawn contour is ~45.3 mjy. Note also that the extended 
galactic emission at 4.9 GHz around SN 2001gd is resolved 
out. Note that there is no contamination of background syn- 
chrotron emission around the position of SN 2001gd, which 
is clearly detected at 4.9 GHz. We note that Nl, SW1, and 
SW2 are also detected at this frequency, despite the 10 years 
elapsed between the 1.4 and 4.9 GHz VLA observations. 

The extension of the 8.4 GHz VLA-D emission from the 



nuclear regions of NGC 5033 (panel c) is about 110" x 50" (~ 
6.9 kpc x 3.2 kpc), a size similar to that of the 4.9 GHz 
radio emitting region. The total flux density for the nu- 
clear region of NGC 5033 above the lowest drawn contour 
is ~1 1.7 mjy. Note that component Nl, as well as the su- 
pernova SN 2001gd, are clearly detected outside the main 
bulk of the nuclear radio emitting region. Note also that a 
new component, El, is detected close to the inner regions of 
the disk, while components SW1 and SW2 are not detected 
above a limiting flux density of 57/iJy/b 

Panels (d) and (e) of Fig. [T] correspond to radio images 
of the nuclear regions of NGC 5033, obtained at 1.5 and 
4.9 GHz from archival observations made with the VLA in 
A- and B-configuration, respectively. The peaks of bright- 
ness are of 5.7 and 3.4 mjy beam -1 at 1.5 and 4.9 GHz, 
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Table 1. VLA sources within the galaxy NGC 5033 



Source a(J2000.0) 5(J2000.0) S[ A S{ 4 Sf g 5| g S£ 4 5| 4 
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* The coordinates correspond to the putative nucleus of NGC 5033, dubbed NGC 5033*, as seen with the VLA-A and VLA-B at 
1.4 GHz (October 1999) and 4.9 GHz (August 199), respectively. and Si are the peak and integrated flux densities for each given 
component, at a frequency v. * Peak and integrated flux densities for NGC 5033, as seen with the VLA-C a nd VLA-D at 1.4 (A ugust 
1993), 4.9 (April 2003), and 8.4 GHz (April 2 003). a Sour ce tentatively identified as background source by iHo fc Ulvestadl J2001T) ; b 
Source position coincident with an H II region (Evans 1996). 



respectively, and correspond to the emission from the core. 
A jet pointing eastwards from the core is also seen at both 
frequencies. The different position angle of the jet, as ob- 
tained at 1.4 GHz with respect to 5.0 GHz cannot be fully 
accounted for by the different synthesized beam at each 
of the two frequencies. In fact, while the 1.4 GHz VLA 
data displays a very well collimated jet, the 4.9 GHz VLA 
data shows a hint for a double-layer jet structure, which 
is very suggestive of the spine/shear layer model proposed 
by lLaing fc Bridl (2002) to explain the synchrotron radio 
emission of Fanaroff-Riley radio galaxies. We also note that, 
while the position angle of the upper layer of the 4.9 GHz jet 
coincides well with that seen at 1.4 GHz, the bottom layer of 
the 4.9 GHz jet is not seen at 1.4 GHz. This could point to 
opacity effects at the base of the jet playing a relevant role, 
which could be responsible for shaping the position angles 
seen at 1.4 and 4.9 GHz. 



3.1 Background sources 

We obtained the 1.4 GHz VLA-C image in Fig. 1 (panel a) by 
using wide-field techniques, to accurately take into account 
the contribution to the total radio emission from confus- 
ing sources within an area of 50' x 50' around the nucleus of 
NGC 5033. We used the AIPS task SAD to search for point- 
like sources above a flux density level cutoff of five times the 
off-source (lrms=93 /iJy beam -1 ). This resulted in the de- 
tection of a total of 49 background sources (i.e., outside the 
lowest drawn contour in panel (a) of Fig. [1} , detected above 
five times the off-source rms (1 rms = 93 /iJy beam" 1 ). In 
Table [2j we list the coordinates, peak, and integrated flux 
densities for each source component detected at 1.4 GHz. 
A number of sources exist in the literature (see caption of 
Table [2]), but many of them are detected for the first time. 
The most strikin g finding among t hese background sources 
is B2 1312+36 jColla et all I1973T ), which showed a peak 
flux density of 226.1 mjy beam _1 on 7 August 1993, and of 
4.6 mjy beam _1 at 4.9 GHz on 8 April 2003. The enormous 
difference in flux density between the two frequencies would 
imply an ultra steep spectrum for B2 1312+36. However, 
the non-simultaneity of these observations do not allow for 
a precise estimate of a. Fortunately, we were able to make 
use of snapshot 1.4 GHz VLA observations taken also on 
8 April 2003, to precisely determine a. The 1.4 GHz peak 



flux density of B2 1312+36 was of 190.3 mjy beam , in- 
dicating otxA = -3.1 ± 0.1 (S v oc i/"), and thus confirming 
that B2 1312+36 is an ultra steep spectrum (USS) source. 
Such steep spectrum is found either in pul s ars, or in very 
varia bles sources (e.g. lDe Breuck et alll2000l llzvekova et all 
Il98lf) . A FIRST image of B2 1312+36, taken on 17 July 
1994, shows a 1.4 GHz peak flux density of 438.6 mjy/b, 
which confirms this source displays strong flux density vari- 
ations, and hence is very unlikely to be a pulsar. The source 
has no optical counterpart at the limit of the POSS I plates 
(m,R «20.0) and SDSS plates (m v «24.0), neither an infra- 
red counterpart at the limit of the 2MASS plates (J~15.8; 
H«15.1; Kf»14.3). 



4 DISCUSSION 

4.1 The radio spectrum and magnetic field of 
NGC 5033 

We used the VLA archival data obtained in 1999 under 
project AU0079 to derive the spectral index of the core- 
jet structure that is seen in panels (d) and (e) of Figure [T] 
Since the angular resolution of the 4.9 GHz VLA obser- 
vations in A-configuration is essentially the same achieved 
by the 1.5 Ghz VLA observations in B-configuration, we 
used the peaks of brightness at each frequency (as given by 
the AIPS task IMSTAT) to estimate the spectral index of 
the compact, nuclear source. We obtained a = 0.44 ± 0.04, 
which is in agreement with the spectral index published by 
IHo fc Ulve stad (2001) and confirms the AGN-like nature of 
the radio emission from the nuclear region of this Sy 1.5 
galaxy. Since the angular resolution of the 4.9 GHz VLA-D 
is very close to that of the 1.4 GHz VLA-C, we can estimate 
the spectral index of the kiloparsec scale disk of NGC 5033. 
Using the peaks of brightness at those frequencies, we ob- 
tain q = —0.93 ± 0.02(S' J / oc u a , indicating an steep syn- 
chrotron spectral index of the radio emission within the in- 
ner few kiloparsecs of NGC 5033. iLisenfeld fc Voelkl [2000 
have shown that galaxies having steep spectral radio in- 
dices (cv < — 0.9) with increasing distance from the disk 
tend to have significant inverse Compton and synchrotron 
energy losses of cosmic ray ele ctrons, indicating that the ir 
escape rate is low (e.g.. M 82 (|Seaquist fc Odegardlll99ll'); 
NGC 253 (|Carilli et all 1 19921 ); NGC 2146 l|Lisenfeld etall 
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Table 2. 1.4 GHz Background sources^ 
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' Background sources of the 1.4 GHz NGC 5033 field, from 
VLA observations obtained on 7 August 1993 with the VLA- 
C. Sf 4 and S[ 4 are the peak and integrated flux densities of 
each component, respectively. We used the AIPS task SAD to 
search for point-like sources above a flux densi ty level cutoff of fiv e 
times the off-source (1 rms= 93 fiJy beam" 1 ). HColla et al.Nl973h: 
^McMahon et al.l J2002t) ; 3Ho fe Ulvestadl <200ll) ; ^Hales et al.l 
il988|); 



Il996t ): NGC 4666 |Sukumar et alJll988h V While the spec- 
tral index we find for the kiloparsec region of NGC 5033 
hints to the fact that energy losses are important, we can- 
not exclude that the observed spectrum is due to a rather 
normal, relativistic electron energy distribution with p w 2.8 
(N(E) oc E~ p ). 

The 1.4 GHz monochromatic luminosity of the core of 
NGC 5033, as obtained with the VLA in A-configuration, is 
equal to of Li.s = (1.2 ±0.2) x 10 27 Df 3 ergs s" 1 Hz" 1 , and 
the integrated isotropic radio luminosity between 300 MHz 
and 30 GHz is, for a = -0.44, of L corc = (1.5 ± 
0.2) x 10 37 £> 13 ergs s . From our VLA observations in D- 
configuration, we obtain a 1.4 GHz monochromatic lumi- 
nosity for the nuclear regions of the galaxy NGC 5033 of 
Li.4 = (3.0 ± 0.3) x 10 28 Z> 2 3 ergs s" 1 Hz" 1 . If the spectral 
index of a = —0.9 extends at least from 300 MHz to 30 GHz, 
then the integrated isotropic radio luminosity for the nuclear 
regions of the galaxy is Lr = (2.3 ±0.2) x 10 38 Dj a ergs s" 1 . 
We can then estimate the average magn etic field in the radi o 
emitting regions NGC 5033. Following fPacholczykl (|l970h . 
we have 

B min = (4.5 c 12 /^) 2/7 (1 + ^) 2/7 i?" 6/7 L R /7 (1) 

where Lr is the radio luminosity of the source, in ergs s" 1 ; 
7? is the characteristic linear size of the source, i n cm; C12 is 
a slow ly- varying function of the spectral index l|Pacholczvkl 
Il970h ; <j> is the filling factor of fields and particles, i.e. the 
ratio of the volume filled with relativistic particles and mag- 
netic fields to the total volume occupied by the source; and 
ip is the ratio of the heavy particle energy to the electron 
energy. Since the value of <j> is highly uncertain, we will set 
it to 0=0.5 for the sake of simplicity. The value of xj) de- 
pends on the mechanism that generates the relativistic elec- 
trons, and can range from 1 to 2000. We will adopt here 
a value of rj) = 100, which seems appropriate for galaxies 
l|Moffetl li"973h . As the characteristic size of the source, R 
we took an angular region of 30"x30"centered in the nu- 
cleus of NGC 5033, which approximately corresponds to the 
inner region of the galaxy disk. The corresponding approxi- 
mate flux density values at each frequency are of 58, 21, and 
6 mjy at 1.4, 4.9, and 8.4 GHz. With the above values, Eq.[T] 
results in minimum magnetic field values of 23, 22, and 18 
/xG at 1.4, 4.9, and 8.4 GHz, respectively. Therefore, if the 
situation is close to equipartition, the average magnetic field 
in the inner region of the disk is of about 20 piG. 

4.2 Radio/FIR correlation and the star-formation 
rate of NGC 5033 

The existing linear correlation between the total radio con- 
tinuum emission and the far- IR luminosity (Lfir) is well 
known in "normal" galaxies (|Condon|[l992l ). and is gen- 
erally interpreted as being due to the presence of mas- 
sive stars. These massive stars provide relativistic par- 
ticles via supernova explosions and heat the interstellar 
dust, which then re-radiates the energy at FIR wavelengths 
(|Helou et al.lll985l ; ICondonlll992l '). The ratio of infrared to ra- 
dio luminosity is u sually expressed through the g-parameter 
jHelou et al.ll 19851 ). which is defined as q = log [(FIR/3.75 x 
10 12 )/Si.4 GHz], where Si.4GHz is the observed 1.4 GHz flux 
density in units of W m~ 2 Hz -1 , and FIR is given by 
FIR = 1.26 x 10" 14 (2.58560 ^m + Sioo Mm) where S 60Alm and 
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Sioo/jm are IRAS 60 /jm and 100 fim band flux densities, in 
Jy (Helou et al. 1988). 

Therefore, q is a measure of the logarithmic FIR/radio 
flux-density ratio, and is an indicator of the relative im- 
portance of an AGN, or starburst in a galaxy. Most galax- 
ies in the IRAS Bright Galaxy Sample have q « 2.34, 
though some galaxies have smaller q values due to addi- 
tional cont ributions from compact radio cores and radio 
jets/lobes (|Sanders fc Mirabell fl996T) . The bulk of the far 
infra-red emission comes from a region of angular size of 
radius 30" (~ 1.9 kpc), as measured from an archival ISO 
image of NGC 5033. Hence, to compute q, we estimated 
the 1.4 GHz radio emission from a circumnuclear region of 
the same angular size, using the AIPS task JMFIT. The re- 
sulting value is of about 58 mjy, which translates into an 
isotropic luminosity of L 1A = 1.2 x 10 21 W Hz' 1 . The IR 
flux densities at 60 and 100 /J,m are Sep um = 16.20 J y and 
5*100 Mm = 50.23 Jy, respectively l|Sanders et al.l 12003). The 
resulting q-parameter -after subtraction of the radio emis - 
sion contributed by the nucleus- is q=2.83. lYun et all (|200ll ) 
pointed out that values of q less than 1.64 would indicate 
the radio dominance of an AGN. Thus, it appears that the 
radio nucleus of the Sy 1.5 galaxy NGC 5033 is of very low- 
luminosity, and the radio emission from its circumnuclear 
region is dominated by starburst activity. We can charac- 
ter ize such a sta r burst by, e.g., following the prescriptions 
by IScoville et al. 

I (|l997h . which characterize a starburst by 
its total luminosity, L, Ly continuum production, Q, and 
accumulated stellar mass, M*. Those values are obtained as 
power-law approximations of the lower and upper mass cut- 
offs for stars, mi and m u , the constant rate of star formation, 
M, and the burst timescale, is. In particular, a starburst 
lasting 10 Myr and producing stars at a rate of 2.8 Mq yr -1 
in the range mi = 1 Mq and m u — 40 Mq , yields essentially 
the same observed far-infrared luminosity of NGC 5033 
{Lfjr — 1.2 x 10 10 Lq), and results into a Ly continuum 
production of Q ta 3 x 10 53 s"\ and M* » 9 x 1O 7 M . To 
estimate the supernova rate, we used a minimum mass for 
yielding Type II supernovae of 8Mq . The corresponding su- 
pernova rate was of 0.045 yr~ 1 , which corresponds to about 1 
SN event every 22 yr. Therefore, we find that a rather mod- 
est (both in time and intensity) starburst scenario is able 
to satisfactorily account for both the observed far infra-red 
and radio luminosity of NGC 5033. (Interestingly, NGC 5033 
has been host galaxy to at least three supernovae in the last 
60 yr: SN 1950C, SN 1985L, and SN 2001gd, although the 
supernova rate we have inferred applies for the nuclear and 
circumnuclear regions, where the main bulk of radio and 
FIR emission come.) 

4.3 The nature of the non-nuclear radio sources 
in NGC 5033 

We now discuss the nature of the four sources (Nl, SW1, 
SW2, and El; see Table [1} detected in the VLA images 
well within the radio contours of the 1.4 GHz VLA-C radio 
emission (panels (a) and (c) in Fig. fJJ 

4.3.1 Source Nl 

Nl is catalogued in lHo fc Ulvestadl (|200ll ) as a background 
source. Nl is also detected in our 4.9 and 8.4 GHz images, 



taken 10 years later than the 1.4 GHz in Fig. 1. (see also 
Table [TJ. Its spectral index between 4.9 and 8.4 GHz is 
a ~ —1-1, and if we assume no variability between 1993 
and 2003, the integrated spectral index from 1.4 to 8.4 
GHz is a = -0.8 ± 0.1. From our 8.4 GHz VLA-D image, 
which yields the highest resolution, the coordinates of Nl are 
a=13 h 13 m 33 , !44±0 ? :02, <5=36°37'28"7±0"2 (J2000.0), These 
coordinates are coi ncident with those reported for FIRST 
J131333.4+363728 (|Becker et al.lll995h , which suggests Nl 
being a background source. Moreover, Nl is about 6" away 
from the closest H II regions, so it is likely to be unrelated 
to them. We therefore conclude that Nl is a background 
source. 



4.3.2 Source SW1 

SW1 had a 1.4 GHz peak flux density of 4.5 mjy on April 
1993, and was also detected at 4.9 GHz and 8.4 GHz (peaks 
of 0.66 mjy and 0.36 mjy, respectively). The source was ten- 
tative ly identified as a background source bv lHo fc Ulvestadl 
|200ll ). at a flux density level of 0.3 mjy at 4.9 GHz. The 
source therefore shows some variability, but the up-and- 
down behaviour at 4.9 GHz is at odds with it being a radio 
supernova. The clo sest H II re gion is about 5"away from 
the radio position ((Evans 1996), which suggests it is unre- 
lated to it. Finally, the source was detected by ROSAT, but 
there is no detection in the optical. We therefore suggest, as 
previous authors, that SW1 is also a background source. 

4.3.3 Source El 

El, detected in our 8.4 GHz VLA-D observations 
(panel (c) in Figure [T|, also coincides with the s ource 
J131330. 1+363537 detected by iHo fc Ulvestadl l|200ll ) with 
a flux density of S*4.g ~ 0.20 ± 0.06 mjy, and which these 
authors tentatively identified as a background source. Our 
continuum 8.4 GHz radio observations indicate an angu- 
lar size for El of ~ 7.4" (corresponding to a linear size of 
~ 470 pc). The 8.4 GHz flux density of above 2 a is Ssa = 
(123 ± 19) fi Jy. We also used XMM-Newton archival data of 
NGC 5033 taken on 18 December 2002 as part of an observ- 
ing program to ob s erve S N 2001gd, and refer the reader to 
IPerez- Torres et all ([2005) for further technical details. We 
analyzed the X-ray data in a region ofz 9"in radius (ap- 
proximately the same size as the radio emitting region) , and 
centered at the peak of 8.4 GHz radio emission. The X-ray 
spectrum has 200±20 counts, and is best fit by an optically 
thin thermal plasma with temperature kT — 2.3lJ'g keV 
and column density Nu = (1.8lJ'g) x 10 21 cm -2 , for a galac- 
tic column density of Ng* = 1.1 x 10 20 cm" 2 . 

For the reasons outlined in the previous paragraph, we 
suggest that El is a radio "spur" which emanates from 
the disk of NGC 5033. Such radio spurs can be observed 
as a consequence of the formation of superbubbles of size 
~ 300 — 1000 pc, if sequential Type II supernovae explo- 
sions occur every ~ 10 5 yr during the lifetime of a giant 
molecular cloud, canonic ally taken to be ~ 1 — 2 x 10 7 yr 
|Norman fc Ikeuchilll989l and references ther ein) and have 
been previously detected, e.g., in NGC 253 (|Heesen et al.l 
2004). The energy input by stellar winds and multiple su- 
pernovae resulting from OB associations results in an over- 
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pressurized region (the bubble), that expands and pops out 
of the disk. The supernovae explosions would then acceler- 
ate the relativistic electrons that, immersed in a magnetic 
field, would be responsible for the observed synchrotron ra- 
dio emission, which is observed as a radio spur. 

The non-thermal spectral index for the spur, a\% ~-0.9, 
is similar to the non-thermal spectra observed for radio su- 
pernovae. If the synchrotron spectrum holds at least between 
300 MHz and 30 GHz, the total radio luminosity of the spur 
would be 1.1 x 10 36 D\ 3 ergs s _1 , a typical value for relatively 
old radio supernovae, or young supernova remnants. We can 
then estimate the minimum total energy in relativistic parti- 
cles and fields, and the approximate equipartition magnetic 
field ( Eq.[l]) in the radio spur from minimum energy argu- 
ments (Pach olczvklll970h : 



l-\ i ; J.3/7 D 9/7 r 4/7 
-Emin = C13 (1 + 1p) ' (f> ' R ' £ R 



(2) 



where C13 is a slowly-varying function of the spectral index, 
and the rest of the parameters have the same meaning as in 
Eq. [T] and we will also use <j> = 0.5 and if) = 100 for the sake 
of simplicity. 

We then obtain -E m m ~ 4 x 10 52 ergs and f? m i n ^ 
13 /iG. The lifetime of the radio spur can be obtained by 
assuming that its observed radio luminosity has been con- 
stant. In that case, r spur = _E rc i/LR,. Since the relativistic 
particle energy is E vc i = 4/7(1 + ip)Em in ~ 4.3 x 10 50 ergs, 
we get T spur ~ 16.1 Myr. This value agrees well with the ex- 
pected lifetime for giant molecular clouds, estimated to be 
~ 10 - 20 Myr. 

The characteristic lifetime of ele ctrons undergoing ra- 
diative synchrotron losses is (see, e.g.. |Pacholczvklll970l ) 



E 



33.4 B 



-3/2 -1/2 



Myr 



(3) 



* yu ~(dE/dt) syn 

where Bio is the magnetic field in units of 10/iG, and v\ 
is the critical frequency, in GHz. Hence, the non-thermal 
electrons emitted in the radio spur have lifetimes of about 
19.0, 8.6, and 6.6 Myr at 1.4, 4.9, and 8.4 GHz, respectively. 
It then follows that while the low-frequency radio emission 
(1.4 GHz) of the radio spur may be due to synchrotron losses 
of the initially injected electron population, the higher fre- 
quency radio emission requires either subsequent injection of 
fresh electrons, or a reacceleration of the old ones, or both. 
These facts are all in agreement with the cosmic ray elec- 
trons in the radio spur being injected by young supernovae 
and supernovae remnants occurring during the lifetime of 
the host giant molecular cloud. 

In summary, while we cannot completely exclude that 
Nl is a background radio source, there is huge evidence for it 
being a galactic source, namely a radio spur. Indeed, the size 
(~ 470 pc), radio spectrum (a ~ —0.9), high temperature 
of the ionized plasma (T c ~ 2.5 x 10 7 K), and radiative 
lifetime of the putative spur (~ 16 Myr) suggest that it could 
be a hot superbubble formed as a consequence of sequential 
supernova explosions occurring during the lifetime of a giant 
molecular cloud. 

4.3.4 Source SW2 

SW2 was detected at 1.4 GHz in August 1993 (S[ A = 
1.85 mjy beam -1 ), and further identified at 4.9 GHz in 
April 2003 (S£ 9 = 0.32 mjy beam" 1 ). The radio location 



of SW2 is spatially coincident , within less than 1.5", with 
an H II region of NGC 5033 l|Evans!ll996T ) and an optical 
coun terpart from the NOMAD catalogue llZacharia s et al.l 
2004! ) (observed in 1963 with B = 19.64 and R = 19.97). In 
turn, this H II region is close to a region detected with the 
WFPC2 camera onboard HST, which suggests their emis- 
sion is physically related. Thus, SW2 would appear to be of 
galactic origin. A possible scenario could be that proposed 
for El in the previous paragraphs, i.e., that SW2 is a ra- 
dio spur. The angular size of the spur, about 8"in radius 
(~500 pc) is a typical value for the putative superbubble. 
An alternative (galactic) scenario could be that of a radio su- 
pernova that exploded around 1993. The observed 1.4 GHz 
luminosity would then correspond to an event about three 
times fainter than SN 2001gd at its peak. However, both the 
radio spur and the radio supernova scenarios are very dif- 
ficult to reconcile with the fact that we did not detect any 
8.4 GHz emission on April 2003 above 57/iJy (3 a) around 
the position of SW2, which implies a very steep spectrum 
a t'.t ~ —3.0, and which rather suggests an extragalactic ori- 
gin for this source. In sum, while it would seem that SW2 
belongs to NGC 5033, we cannot rule out the possibility that 
this is a background source. 



5 SUMMARY 

We have presented continuum VLA observations of the 
Sy 1.5 galaxy NGC 5033, made at 4.9 and 8.4 GHz on 8 
April 2003, and also archival VLA data at 1.4 and 4.9 GHz, 
which probe the radio emission of this galaxy from regions 
of less than a hundred parsecs up to more than 35 kpc in 
size. We summarize our main results as follows: 

• The high-resolution VLA images show a core-jet struc- 
ture for the Sy 1.5 nucleus. The core has a moderately steep 
radio spectrum, (S„ cx u a ; aft = -0.44±0.04). The 1.4 GHz 
radio emission, as traced by the VLA in C configuration, 
correlates very well with the optical emission of the whole 
galaxy, and delineates exquisitely the spiral arms. Our 4.9 
and 8.4 GHz simultaneous observations with the VLA-D 
show only significant radio emission from the inner disk of 
the galaxy. Combining the data at 1.4, 4.9, and 8.4 GHz, 
we find that the radio emission of the inner disk of the 
galaxy shows a steep spectrum (af t ~ —0.9), which is usu- 
ally found in galaxies where radiative losses are high in the 
disk of the galaxy, so that the escape rate of their cosmic 
ray electrons is low (Lis enfeld fc Voe"lk1l2000h . 

• If the synchrotron spectrum extends at least from 
300 MHz to 30 GHz, the isotropic radio luminosity 
of the inner disk of NGC 5033 is L R (2.9 ± 0.3) x 
10 21 D 2 3 Wm" 2 Hz~\ Combining this value with the far- 
infrared luminosity of the galaxy (Lfir = 1.20 x 10 L0), 
we obtained a value for the q parameter of 2.83, which in- 
dicates that the radio emission from the inner regions of 
NGC 5033 is mainly powered by a recent starburst. In fact, 
the compact core-jet structure at the center of the galaxy, 
contributes only about 7% of the total 1.4 GHz radio lu- 
minosity. We find that a relatively short starburst, lasting 
10 Myr and and producing stars at a rate of 2.8M0yr _1 
in the range mi = 1 Mq and m u = 40 Mq , reproduces very 
well the observed radio and far-infrared luminosities, and 
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results into a supernova rate of 0.045 yr -1 . This supernova 
rate corresponds to about 1 SN event every 22 yr. 

• We find evidence for the existence of a radio "spur" (El) 
of radius r ~ 470 pc, located at a distance of about 2.1 kpc 
above the disk plane. Its radio spectrum (a ~ —1.1), high 
temperature of the ionized plasma (T e ~ 2.5 x 10 7 K), and 
lifetime (r ~ 16 Myr) strongly suggest that it is a hot su- 
perbubble formed as a consequence of sequential supernova 
explosions occurring during the lifetime of a giant molecular 
cloud. 
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